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Abstract 

A Ward-Takahashi identity, as a consequence of gauge invariance and in 
a form that relates self-energy to the two-particle Bethe-Salpeter scattering 
kernel, was first derived by Vollhardt and Wolfle for a system of independent 
particles moving in a random medium. This is generalized to a class of in- 
teracting electronic systems in materials with or without random impurities, 
following a procedure previously used for classical waves transport in disor- 
dered media. This class of systems also possesses other symmetry properties 
such as invariance under time translations and local spin rotations, which im- 
ply local conservation laws for energy and spin current. They imply additional 
Vollhardt- Wolfle type identities. We present non-perturbative derivations of 
these identities, and consider the constraints they impose on the relationship 
between the self-energy and the two-particle scattering kernel. 
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Ward identities are statements of local conservation laws, and reflect the corresponding 
symmetry properties of the system. The original Ward identity [T|] and its generalization by 
Takahashi || were derived for quantum electrodynamics, as a consequence of local current 
conservation, and relate the current vertex function to the electron propagator. In investigat- 
ing the problem of the Anderson localization, Vollhardt and Wolfle H considered a system 
of independent electrons moving in a random potential, and derived a Ward- Takahashi iden- 
tity in a form that relates the electron self-energy to a two-particle Bethe-Salpeter kernel 
induced by statistical randomness of the medium. More recently, similar identities were 
derived for classical waves transport in disordered media |4[], for which, however, energy 
conservation underlies the physics involved. It is clear from the latter derivations, which 
make use of a Bethe-Salpeter type integral equation, that Vollhardt and Wolfle's version 
of the Ward identity is valid for a considerably wider class of systems, such as interacting 
electrons moving in materials with or without random impurities. For this general class 
of electronic systems, energy and spin, in addition to charge, are also conserved, reflecting 
invariance properties of the system under time translation and spin rotation, and implying 
additional Vollhardt-Wolfle type identities. We consider in this paper the derivations of 
these identities, and their implications on the self-energy. 

The model Hamiltonian we consider is of the form, 



H = I d 3 x 



— V *(x, t)- V $(x, t) + Vi(x)^(x, t)^(x, t) 
2m 

+ d 3 xd 3 x'ty(x, t)^f(x, t)V 2 {x - x')^(x ,t)m(x' : t), 



(1) 



where \l/ is the two-component electron field, quantized according to the usual equal-time 
anti-commutation relations, V\ the potential provided by atoms in the material, including 
randomly distributed static impurities, and V 2 the interaction potential between any two 
electrons with the property 

V 2 (x-x , ) = V 2 (x' -x). (2) 
The corresponding Lagrangian density is given by 



L = ity(x) * (x) - ^— V *(a;)- V - V^(x)*(x)*(x) (3) 



For convenience, we use the four-dimensional notations with a; M = (t,x,y,z) and = 
(1,-1,— 1,-1) for // = 0,1,2,3. Due to the presence of the non-local interaction term 
V2 in the Lagrangian density, one should exercise care in deriving the equation of motion. 
Instead of naively applying the Euler-Lagrange equation, one considers the minimization of 
the action, S J Ld 4 x, to derive the expected equation of motion : 

i (~^V 2 + k) * + J d 3 x'y(x',t)y(x',t)V 2 (x-x')y(x). (4) 

The Lagrangian density is invariant under the gauge transformation, 

— ► e iX ^(x), (5) 
— ► e~ iX ^(x), 

the time translation, 

* — >l/ = t + e, (6) 
— ► V(x,t' -e), 
V(x,t) — >^(x,t' -e), 

and the spin rotation 

— ► e <3 -*tf(a;), (7) 

where <r are the Pauli spin matrices. The corresponding conserved quantities are charge, 
energy, and spin, respectively. 

In disordered media, there is a random distribution of impurities. We need to consider a 
statistical ensemble of impurity configurations and average over the ensemble. Under rather 
general assumptions on the statistical ensemble, as has been fully discussed by, e.g., Frisch 



||, statistical averaging results in an effective electron-electron interaction, which adds to 
the real electron-electron interaction potential V2. 

Local current conservation as a consequence of gauge invariance is of the form 

d,f = 0, (8) 

where 

f = j k = ^-[m k q - (d k y)^}. (9) 

It implies that 

d?{0\T{f(x)*(y)V(z)}\0) (10) 

= -5\x - y)(0\T{y(x)V(z)}\0) + 5 4 (x - z){0\T{V(y)V(x)}\0) 

= -S 4 (x - y)-G(x, z) + 5\x - z)-G(y, x). 
1 1 

The above equation holds for any single configuration of the impurities. When statistical 
averaging on an ensemble of impurity configurations is applied to it, an additional statistical 
interaction infuses itself into the structure of the expressions on both sides of the equation. 
An irreducible current vertex function is defined by 

(0\T{ 3fl (x)^(yW(z)}\0r = J d%d'z X -G^\y - V )T ^xlQ-^ (C - z), (11) 

where expressions on both sides of the equation refer to configurationally averaged quantities, 
and use has been made of the fact that averaging restores translational invariance and, as a 
consequence, the averaged propagation function, G^ e \ depends only on the difference of the 
two coordinates. Averaging applied to both sides of (10) leads to 

an f rfV 4 C-c (e) (?/ - v)rM*\0-G {e \( - z) (12) 

= -5\x - y)-G {e \x -z) + 6\x - z)-G {e \y - x), 
i i 

which implies the Ward-Takahashi identity : 

d^T^r)\x\C) = G^-\r] - x)~5 4 {x - () - S 4 (r] - x)\g^~\x - Q. (13) 
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For a system of non-interacting electrons moving in disordered media, Vollhardt and 
Wolfle ||, using diagram considerations, derived their version of the Ward-Takahashi iden- 
tity that relates the electron self-energy with the two-particle Bethe-Salpeter kernel instead 
of the current vertex function. Rather than following their analysis in terms of diagrams, 
we here adopt the non-perturbative procedure previously used for the derivation of energy- 
conservation Ward-Takahashi identities for classical waves transport in disordered media 

I- 

A crucial input in the non-perturbative derivation is the following integral equation for 
the vertex function T M , 

r^faixlo-rjWfalilc) (14) 

= / d'ydyd'zd'z'-G^iy'-y^^ylx^-G^iz-z^K^C^'z'), 

where K is the Bethe-Salpeter kernel, and denotes the vertex function corresponding to 
a reference homogeneous medium and in the absence of the interaction V2. It is of identical 
form as the one in quantum electrodynamics ||. The validity of this integral equation, 
schematically illustrated in Fig. 1, is on the same footing as the Bethe-Salpeter equation. 
We note that all quantities in the above integral equation refer to statistically averaged ones. 
Combining the two equations (13) and (14) results in the Ward identity in the Vollhardt- 
Wolfle form : 

£(77 - x)5\x - C) - S(ar - ()5\v ~ x) (15) 
#yd*z[8\y - x)G ie \x -z)- 5\x - z)G {e) {y - x)]U(z -Q^-y^y-z), 



where the self-energy S is defined by 

Q(e)-1 = G (0)-l _ S (16) 

and U is an alternative expression for K, with specified coordinate dependence. The depen- 
dence of K on the coordinates is inferred from the following two properties. First, because 
of the restored translational invariance due to averaging, only three of the four coordinates 



in K are independent, and are chosen to be z — (, rj — y, and y — z. Second, since the two 
channels are mediated by the effective statistical interaction and V2, both of which are static, 
(y — z) can actually be replaced by (y — z). In momentum space, the identity takes the form 

V(q + p)-H(q) = J -^- 3 [G (e) (q + <3) - G (e) (« + P + ")P(q, Q + P, -$)■ (17) 

This is the Ward identity derived by Vollhardt and Wolfle [Q, but here extended to include 
the non-local interaction V2 between any two electrons. 

The Lagrangian density is invariant under the time translation. This implies the local 
conservation law 

d^o = 0, (18) 



where 



T oo(x) = 7T~ V *(x)- V *(a:) + ^i(f)*(x)*(x) 
2m 

+ d 3 x^(tf ) t)^(x',t)V 2 (x-x')y(x)y(x), (19) 
T k0 (x) = ^[d k ~*(x)d y(x) + da*(x)d k y(x)]. 

In particular, the energy density is Too- Similar to the case of current conservation, we con- 
sider in the present case the divergence of the three-point function (0\T{T fl0 (x)^(y)^(z)}\0). 
To proceed, we need to know the equal-time commutation relations [T 00 (x), ^ (?/)], etc.. Due 
to the presence of the non-local interaction V2, we consider the commutation relation with 
the integrated total energy: 

d 3 xT 00 (x,t),V(y,t)] (20) 
J d 3 x5 3 (x-y) I -^ v2 + y i( f ) * + J <Px'^(x" ,t)^{x\t)V 2 {x-x)m(x)^ 



— I Av - [S 3 (x -y) 
Im J 



2m 

There is some uncertainty as to what to do with the 3-dim divergence term in the above 
equation. We have considered various arguments for or against its inclusion in the local 
equal-time commutation relation. The most convincing argument against its inclusion comes 
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from an exercise in deriving the Ward identities by using the path-integral formalism. We 
thus conclude that, on account of (4), 

[Too(x,t)MM] = S s (x-y) M> (x,t), (21) 
and similarly for Consequently 

%(O\T{T0(x)*(y)9(z)}\O) = -5\x - y)^G{x,z) - 5\x - z)-^G{y,x). (22) 
Defining the energy vertex function according to 

(0|T{T^(s)¥(y)*(z)}|0r = / e*V 4 C-G (e) (y - rfiTf \r 1 \x\Q-G^{C - z) (23) 
we obtain from (22) the relation 

d^fHAQ = G^-\v - - + G®-\x - O^/iv - x). (24) 

which is the Ward-Takahashi identity for energy conservation. Compared with the current- 
conservation Ward-Takahashi identity (13), there are here in (24) extra linear time deriva- 
tives, which are the unmistakable signatures for energy. The two identities describe charge 
flow and energy flow, respectively. 

To derive the corresponding Vollhardt-Wolfle type identity, we again make use of an 
integral equation for the energy vertex function, which is identical in structure to the charge 
counterpart, 

rf>fa|s|c)-rjP>fa|s|c) 

= / d 4 yd 4 y'd 4 zd 4 z'-G (e \y' - y)Vf \y\x\z)-G^\z - z>)K( V (; y'z'). (25) 

Combining (24) and (25), we obtain the desired identity 

Efo - x)-^5\x - C) + E(x - O-^-fiv ~ x) (26) 
= J d 4 yd 4 z\b\y - x) w G^\x - z) + 5 4 (x - z^&'Hv ~ x)\U{z -^ v - y ,y- z) 

In momentum space, it is of the form 



q Z(q + p) - (g +Po)S(g) 
= J j^[qoG^(q + Cu)-(q +p )G^\q + p + Cu)]U(q,q + p,-Co). (27) 

The Lagrangian density (3) is also invariant under local spin rotations, for which 

— >e iu} -H(x), V(x) — >^(x)e~ iu} -i (28) 

This symmetry property implies the local conservation law 

8^ = 0, a =1,2,3, (29) 

where the spin current is given by 

S° a = *y *, (30) 

2m v 2 2 ; 

Using a procedure in complete parallel to the cases of gauge and time-translation symmetries, 
we obtain the following Vollhardt-Wolfle type identity for a = 1, 2, 3, 

Efo - *)y * 4 (* - C) - * 4 fo - x)^S(x - C) (31) 
= J d 4 yd 4 z[S\y - x)^\x - z) - 5\x - z)G^\y - x)^]U(z -(, v -y,y-z) 

In momentum-space representation, it is of the form 

£(<7 + p)y-y £(<?) = J^[^G^(q + u)-G^(q + p + u)^]U(q,q + p,-u). (32) 

The collection of the three Vollhardt-Wolfle identities (15), (26), and (31), or their 
momentum-space counterparts (17), (27), and (32) set fairly stringent conditions on the 
relationship between the self-energy and the two-particle scattering kernel. Combining (17) 
and (27) yields the surprisingly simple but strong relation 

£(?) = / ^G^(q + p + u)U(q,q + p,-a). (33) 
In configuration space, it is of the form 



5 4 (r, - x)E(77 - C) = / d 4 yd A z8\x - z)G^\y - z)K( V (; yz). (34) 

We have checked and verified this relation in terms of low-order diagrams. According to 
(33) or (34), all self-energy diagrams can be generated from the corresponding irreducible 
kernel diagrams. For example, the self-energy diagram in Fig. 2a can be generated from the 
irreducible kernel diagram in Fig.2b by connecting with a propagator as shown in Fig. 2c. 
In the verification of the relation (34), we have disposed of the tadpole type self-energy 
diagrams, such as in Fig. 3a, which can not be generated according to (34). The disposal of 
these tadpole type diagrams is based on the following considerations. First, in the traditional 
formulation of the Lagrangian field theory, the Lagrangian density is defined in terms of 
normal ordered products of the field operators. All tadpole diagrams, by the very definition 
of the normal ordering, do not appear at all J?J. Second, if for any reason, normal ordering 
of the products is not adopted, as is effectively the case in the path-integral formulation 
of quantum field theory, tadpole diagrams do appear. However, their contribution to the 
self-energy is "eliminated" order by order by the "mass" insertion diagram in Fig. 3b, due to 
the presence of the mass counter term, actually the chemical-potential counter term in the 
present context, which is introduced in the Lagrangian density to achieve renormalization. 
We note that renormalization counter terms are introduced in either of the above two for- 
mulations, although details of their functions differ somewhat. For example, in the case on 
hand, the "mass" counter term absorbs the (infinite) tadpole contribution to the self-energy 
in the path-integral approach, while in the traditional approach of Lagrangian field theory, 
because the tadpole diagrams do not appear, there is no need for the "mass" counter term 
to absorb them. 

The strong relationship (34) perhaps reflects the fact that both the charge and the 
energy are exclusively carried by the same carriers, the electrons. It would be of interest to 
investigate additional implications of the set of three Vollhardt-Wolfle identities discussed 
in this note. 
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FIG. 2. Schematic relationship between £ and K. 
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FIG. 3. "Tadpole" self-energy. 
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